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' Abstract. We study the Kleinberg problem of navigation in Small World networks 

when the underlying lattice is stretched along a preferred direction. Extensive 
£^ \ simulations confirm that maximally efficient navigation is attained when the length 

G ■ r of long-range links is taken from the distribution P(r) ~ r _Q , when the exponent 

c/2 . a is equal to 2, the dimension of the underlying lattice, regardless of the amount of 

anisotropy, but only in the limit of infinite lattice size, L — > oo. For finite size lattices 
we find an optimal a(L) that depends strongly on L. The convergence to a — 2 as 
L — » oo shows interesting power-law dependence on the anisotropy strength. 



PACS numbers: 89.75.Hc 02.50.-r, 05.40.Fb, 05.60.-k, 
Introduction 
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The small-world phenomenon is one of the most intriguing properties of human society. 



This describes the fact that unrelated people in a society, who are a very large 
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geographic distance apart from one another, tend to be connected by surprisingly short 
chains of acquaintances. This phenomenon was hypothesized in 1929 by Hungarian 
author Frigyes Karinthy [TJ [2] and was first observed experimentally in the 1960's with 
sociologist Stanley Milgram's seminal experiments [3] wherein randomly chosen people 
were selected to mail a letter to an unknown target person, but were only allowed to 
send the letter to a friend, who would pass the letter along to another friend, etc., until 
the target was reached. Successful transmissions took surprisingly few intermediate 
people, lending credibility to the turn of phrase 'six degrees of separation,' popularized 
by Karinthy. Understanding this phenomenon is an important sociological problem. 

To study the underlying mechanism that led to Milgram's results, computer 
scientist Jon Kleinberg modeled a society as follows [U [5] . Begin with a large, regular 
square L x L lattice. Each node is connected to its nearest lattice neighbors and to 
a single random node a large distance away. The probability of nodes i and j being 
connected by such a long range contact is 

^(«) = ^7E^ (!) 

where is the Euclidian distance between the two nodes and the sum runs over all nodes 
in the network except i. Physically, the local lattice connections represent associations 
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with immediate neighbors, fellow townspeople, etc., while long-range contacts might 
model friends or relatives in another city or country. 

We seek to pass the message from a random starting node s to a random target t. 
Of great importance is the fact that each node has no information beyond the locations 
of its contacts and the target node t, so the operational algorithm must be local in 
character. Kleinberg has proved that no local algorithm can do better, functionally, than 
the greedy algorithm [I]: Each message holder passes the message along to whichever 
of its contacts is closest to t, until the message reaches the target. Moreover, for a ^ 2 
the delivery time T (number of intermediate steps) scales as a power of L, while small 
world behavior and the weakest dependence on lattice size emerges for a = 2, where 



The distribution of nodes on a regular square lattice is too rigid, failing to mimic 
important features of actual distributions of populations (or computer routers, etc.). 
In an effort to account for these we have studied Kleinberg navigation in fractals [6], 
showing that the optimal long-contact exponent is then a — df, the fractal dimension 
of the lattice. In this letter we study the effects of anisotropy — another commonly 
encountered distortion of the ideal Kleinberg lattice. Our results indicate that in the 
limit of lattice size L — > oo the optimal contact exponent for two-dimensional lattices 
is still a = 2. For finite L, we find an optimal exponent a(L) quite different from 
the infinite limit. The convergence to a = 2 as L —>■ oo shows interesting power- law 
dependence on the anisotropy strength. 

Anisotropic lattices 

We wish to study the isolated effect of anisotropy on Kleinberg navigation. To do 
this, we begin with a regular square lattice (d = 2) and introduce one of two forms of 
anisotropy. 

• Lattice Anisotropy: The underlying lattice is stretched horizontally, along the 
x-axis, by a factor b > 0, such that the area of each cell goes from 1 x 1 to b x 1. 

• Angular Anisotropy: Long-range contacts are chosen preferably along the 
vertical direction, by a factor b > 0. More precisely, the random angle 9 of each 
long-range contact vector (measured counter-clockwise, from the x-axis) is modified 



Both of these types of anisotropy tend to favor connections in the vertical ^-direction, 
if b > 1, and along the x-direction, if < b < 1. 

Simulations 

To simulate Kleinberg navigation efficiently, we use several tricks and approximations. 
First, rather than testing a finite-size square L x L lattice, we consider an infinite 



T ~ In 2 L . 



to 9': 




(2) 
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lattice and place the source and target at distance L from one another. Since the 
message always progresses toward the target, by the greedy algorithm, the message 
holder remains within a disc of radius L centered on the target node, so in practice only 
a finite number of sites would be explored anyways. Second, the computation of the 
normalizing sum in the denominator of Eq. (OQ) is dependent (in finite lattices) upon the 
location of the node i, and can be very time consuming. The infinite lattice circumvents 
this problem, as the normalizing constant is the same for all nodes. Note, however, that 
J2k r ik a does not converge for a < d. In that case, we imagine a lattice larger than the 
L-disc of activity, with periodic boundary conditions, such that the normalizing factor 
is still the same for all sites. 

Because of the monotonic progression toward the target no site is ever revisited in 
the process. Moreover, as observed by Kleinberg [5], one can think of the long-contact 
link out of node i as being created at the very instant that the message arrives at i. Thus, 
the full lattice is unnecessary, and we need keep in memory only the current location of 
the message holder (and the location of the target). When the message arrives at i, we 
create a random long contact, compare the distances of all five neighbors of i (the four 
lattice neighbors and the long contact) to the target, and move the message to the site 
closest to the target. 

The long contact is created by choosing a random distance and angle, (r, #). In 
order to reproduce the correct P{rij) ~ r^ a the distance r is taken from the distribution 
P(r) ~ r _a_1 , to account for the linear growth of the area of the ring where the contact 
might fall. The angle 9 is distributed uniformly between and 2ir. In the case of angular 
anisotropy, 9 is replaced by 9', according to ([2]). Finally, a vector (r,9) is drawn from 
site i and the contact is placed on the site j closest to the vector tip. 

Because of the anisotropy, the angular displacement from the source to the target 
makes a difference. It is sufficient to test only the two extremes of 9 = and 9 = n/2 
where the target is either parallel or perpendicular to the anisotropy direction. We note, 
however, that anisotropy strength b and a target at 9 = is equivalent to anisotropy 
1/6 and target at 9 = ir/2. For this reason, we simply set the source and target at (0, 0) 
and (L, 0), respectively, throughout, and let b vary both below and above the isotropic 
divide of b = 1. 

Simulations were performed for various values of b over a large range of a and L, 
each averaged 1000 times. For each b and L, the minimum a was computed by first 
fitting a fifth-order polynomial to the averaged data, then using Newton's Method on 
the polynomial's derivative. Finally a m i n was plotted as a function of l/ln 2 L for each 
chosen value of b. These are shown in Fig. [T] and indicate that a m i n — > 2 as L — > oo, 
regardless of b. In Fig. [2] we show detailed results of the extrapolation to L — > oo. 

To further clarify the behavior shown in Fig. [TJ the following procedure was 
performed. First fit a cubic polynomial pb, using least squares, to each 6's curve. Then, 

| A parabola could be fitted to the data closest to the minimum, but we must first know what is 
'closest.' A higher-order polynomial overcomes this difficulty, similar to including higher order terms 
in a series expansion near the minmum of a function. 
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Figure 1: Simulations for lattice and angle anisotropy. A horizontal scale of l/m 2 L 
is used throughout. All curves approach a(oo), regardless of 6. For the lattice case, 
there is a crossover effect where curves for b > 1 dip below the 6 = 1 curve. For the 
angular case, curves for b > 1 approach the infinite limit at differing rates, while curves 
for b < 1 evetually collapse onto the 6=1 curve. These phenomena are further explored 
in Fig. [31 See Fig. [2]for the extrapolated a(oo). 



subtract that polynomial from the isotropic case, pb — p\. This maps b = 1 to the 
horizontal axis and gives the behavior of the b ^ 1 curves "relative" to the isotropic 
curve. These are shown in Fig. [31 The different behavior for each type of anisotropy 
is clear: for both types of anisotropy the results for b > 1 show dramatic differences 
from the isotropic case of b = 1 (the differences for b < 1 and large L are negligible). 
For lattice anisotropy, the b > 1 curves start above the 6 = 1 curve and cross below 
until they eventually converge at a similar rate, as L —>■ oo. For the angular anisotropy, 
the 6 > 1 curves approach a(oo) at a different rate than the 6=1 curve, resulting in 
distinctly different slopes in the plots of Fig. [TBI 

The observed "crossover" behavior present in the lattice anisotropy is somewhat 
unexpected. The crossover point, -^crossover (6), is explored by finding the zero of each 
Pb — Pi- These are plotted in Fig. 0k, and seem to indicate a power law relationship, 
-^crossover (6) ~ 6 2 . Likewise, the different slopes for angular anisotropy, plotted in Fig. [4b, 
show power-law behavior and seem to increase roughly as 6 1//4 . What is responsible for 
these phenomena remains an open question. 

Conclusions 

We have shown, by extensive numerical simulations, that Kleinberg navigation in two- 
dimensional lattices with two types of anisotropy displays the same gross characteristics 
as navigation in isotropic lattices. In particular, the optimal long-contact exponent in 
the limit of infinitely distant source and target remains a = 2, even in the presence of 
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Figure 2: Extrapolation results for lattice (a) and angular (b) anisotropy. Extrapolating 
to 1/ In 2 L — ► with a linear least squares fit to the curves in Fig. [1] shows excellent 
convergence of a(oo) to the expected value of d = 2. Good values should occur when 
the curves are flattest, which happens roughly around 0.25. A more robust fitting 
procedure could be used, but the accuracy of these results imply that it is unnecessary. 
The horizontal lines at a = 2 provides a guide for the eye. 
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Figure 3: Curves relative to the isotropic divide b = 1, for the lattice (a) and angular (b) 
case. To provide a measure of smoothing, cubic polynomials pb were fitted to the curves 
in Fig. [H To clarify the impact of anisotropy, we show the behavior relative to the 
isotropic case, by subtracting pi from each pi,. This maps the isotropic curve to a 
horizontal line and introduces only minor distortion. The crossover behavior for b > 1 
is clearly displayed. 



anisotropy. 

It is worthwhile to note that the actual values for the optimal exponent ct(L) for 
finite L can differ considerably from the limit a = 2, even for reasonably large lattices. 
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Figure 4: Dependence of results on anisotropy strength b for lattice (a) and angular (b) 
anisotropy. The straight lines are of slope 2 and 1/4, respectively. 



Thus, for practical applications the optimal exponent ought to be evaluated on a case 
by case basis. 

The modes of convergence to the limit L — > oo show intriguing power-law 
dependence upon the strength of the anisotropy. A theoretical explanation for this 
behavior remains the subject of future work. 

Acknowledgments 

We gratefully acknowledge NSF award PHY0555312 for partial funding (DbA), and the 
support of a NSF Graduate Research Fellowship (JPB). JMC has been supported as a 
summer research undergraduate student by the McNair program at Carkson University. 

References 

[1] A.-L. Barabasi. Linked: How Everything Is Connected to Everything Else and What It Means for 

Business, Science, and Everyday Life. Plume Books, April 2003. 
[2] Frigyes Karinthy. Chains. In Everything is Different. Budapest, 1929. 
[3] S. Milgram. The small-world problem. Psych. Today, 2:60-67, 1967. 
[4] Jon Kleinberg. Navigation in a small world. Nature, 406(6798), August 2000. 
[5] Jon Kleinberg. The small-world phenomenon: An algorithmic perspective. In Proceedings of the 

32nd ACM Symposium on Theory of Computing, 2000. 
[6] Mickey R. Roberson and Daniel ben-Avraham. Kleinberg navigation in fractal small-world 

networks. Phys. Rev. E, 74(1):017101, 2006. 



